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Abstract 

Neutrino oscillations are examined under the broad requirements of Poincare-invariant scattering 
theory in an ^-matrix formulation. This approach can be consistently applied to theories with 
either field or particle degrees of freedom. The goal of this paper is to use this general framework 
to identify all of the unique physical properties of this problem that lead to a simple oscillation 
formula. We discuss what is in principle observable, and how many factors that are important in 
principle end up being negligible in practice. 
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I. INTRODUCTION 



There is now strong experimental evidence that neutrino flavor eigenstates are mixed via 
a non-diagonal matrix that connects them to neutrino mass eigenstates [1, 2]. The fact that 
space-time propagation is governed by mass rather than flavor eigenstates gives rise to the 
possibility of oscillation between flavor states over space and time. 

There is a large body of literature addressing various aspects of the quantum mechanics 
of neutrino oscillation. Many of the primary issues were set forth by Kayser [3] in 1981. 
Since then, this subject has been examined from a variety of perspectives [4-13] . 

The phenemenon of neutrino oscillations seems simple enough from the perspective of 
introductory quantum mechanics, for which there are countless examples of oscillations in 
two- and three-level systems. There are, however, differences in the neutrino case that tend 
to work against the intuition of the standard examples: the initial conditions for neutrino 
production in reactors or the Sun are generally not controlled; the coherence length over 
which states are virtual can be kilometers rather than subatomic distances, and this affects 
considerations such as what is large or small in a calculation; one of the weak interactions 
that determine the oscillation interval takes place inside the volume of the detector rather 
than far from it; the neutrino kinematics are ultrarelativistic, but not fully so. There is 
also the distinction between what is measured in an experiment and features of a specific 
theoretical approach. These issues have led to a variety of approaches as cited above. Some 
of these approaches have differing perspectives [10, 11], yet all lead to the same simple 
oscillation formula. 

In this paper, we provide an approach to neutrino oscillations from the perspective of 
scattering theory in relativistic quantum mechanics. The virtue of this approach is its 
generality and its focus on what can actually be measured in oscillation experiments. The 
scattering theory approach advocated in this paper has been examined [14] using quantum 
field theory, but the conclusions that follow from our approach are not limited to a field 
theoretic treatment. There are many papers on neutrino oscillations, and many of them 
have significant overlap with this work. Our goal here is to provide a very general framework 
utilizing the S matrix. 

The S matrix for a neutrino oscillation experiment is the probability amplitude for transi- 
tions between states where both initial and final particles are localized near two macroscop- 
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ically separated space-time points; the points of neutrino creation and neutrino absorption. 
We discuss the additional physical conditions that must be satisfied for the transition ampli- 
tude which, when squared, leads to a neutrino oscillation formula. There is also a freedom 
in the choice of variables used to label the single-particle intermediate states, the choice 
of scattering equivalent Hamiltonian, as well as the choice of kinematic symmetries of the 
interaction. Different choices lead to the same scattering matrix, but will normally lead to 
different formulae in specific implementations and different interpretations of the dynamics. 
We find that the conditions that lead to a simple oscillation formula also combine to ren- 
der as insignificant many factors that might otherwise serve to distinguish among various 
theoretical approaches. 

II. FORMAL BASIS: POINCARE INVARIANT QUANTUM MECHANICS 

We present here an overview of the basic ingredients of quantum mechanical systems that 
satisfy relativistic invariance. Further details can be found in Ref. [15]. 

Relativistic invariance in a quantum theory means that a change of inertial coordinate 
system is a symmetry of the theory. In special relativity, inertial coordinate systems are 
related by space-time translations and proper orthochronous Lorentz transformations. The 
group generated by these transformations is the Poincarc group. Proper orthochronous 
Lorentz transformations do not include the discrete transformations associated space refiec- 
tion and time reversal, which are broken by the weak interaction. Wigner proved [16] that 
the existence of a unitary representation, U (A, a) of the Poincare group is both a necessary 
and a sufficient condition for a quantum theory to be relativistically invariant. 

The dynamics are given by the time-translation subgroup of the Poincarc group, which 
is generated by the Hamiltonian, H. Consistency of the initial value problem requires that 
interactions must also appear in additional Poincare generators. This means that at most a 
subgroup of the Poincare group can be independent of interactions. When such a subgroup 
exists, it is called a kinematic subgroup. While the existence of a kinematic subgroup is 
a choice of representation, the choice of representation may imply kinematic symmetries of 
the interaction that have no impact on the scattering matrix elements; this is an example of 
a representation-dependent feature of an interaction that is not experimentally observable. 

The largest kinematic subgroups were classified by Dirac [17]; the classification was com- 
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pleted by Patera and Winternitz [18]. The three largest kinematic subgroups are the three- 
dimensional Euclidean group (instant-form dynamics), the Lorentz group (point-form dy- 
namics) and the subgroup that leaves a three-dimensional hyper-plane tangent to the light 
cone invariant (null-plane dynamics). 

Interactions in a dynamical model where the kinematic subgroup is the three-dimensional 
Euclidean group are translationally invariant and thus conserve three-momentum; interac- 
tions in a dynamical model with a null-plane kinematic symmetry are invariant with respect 
to translations in the null plane and thus conserve components of the four momentum that 
the generate translations on the null plane; while interactions in a dynamical model where 
the kinematic symmetry is the Lorentz group are Lorentz invariant. 

Scattering theory in a Poincare invariant quantum theory can be formulated using the 
same time-dependent methods, based on dynamical and asymptotic Hamiltonians, H and 
Hq. that are used in both non-relativistic quantum mechanics and quantum field theory. 

Scattering states \^{t)) are solutions of the time-evolution equation (Schrodinger equa- 
tion) 

|^±(t)) = e-^^*|^±(0)), (1) 

where the initial condition, |^±(0)), is determined by an asymptotic condition. There are 
two natural scattering asymptotic conditions; they require that the state |^±(i)) approach 
a state of non-interacting particles, |$±(i)), in the asymptotic future or past: 

hm |||*±(t))-|$±(t))||=0 (2) 

t—>-±oo 

where 

|^±(i)) = e-^°*|$±(0)) (3) 

is the non-interacting state. In the field theoretic case the Hq must include self-interactions. 
The appropriate generalization of the asymptotic condition (2), first formulated by Haag 
and Ruelle, is discussed in [19] [20] [21]. 

Equations (1-3) imply that the initial condition for the scattering states are related to 
the initial conditions for the non-interacting states by 

|*±(0))= lim e^^*e-^-^°*|$±(0)) :=Q±|$±(0)). (4) 

t— >±oo 

The probability amplitude that a state prepared to become |$-(i)) in the asymptotic 
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past scatters into a state that becomes in the asymptotic future is 

{S) = {^+m-{t)) = (^+(0)|^_(0)) = ($+(0)|5|$_(0)) (5) 

where 

S := (6) 

is the scattering operator, and we have exploited the invariance of the probabihty amphtudes 
under time translation to emphasize that the probability amplitude can be computed using 
states at any common time. We will use this in formulating scattering involving neutrino 
intermediate states. 

We note that if H is transformed with a unitary transformation A such that H' — AHA' , 
and A satisfies 

hm ||(/-^)e-^°*|*)|| = (7) 

t—^±oo 

for both time limits then 

S = S{H,Ho) = S{H',Ho). (8) 

Note that Hq is not transformed, so the representation of the asymptotically free particles 
remains unchanged. This means V — H — Hq and V — H' — Hq are distinct interactions that 
cannot be experimentally distinguished. Operators A satisfying the asymptotic property (7) 
are called scattering equivalences. They preserve the scattering matrix without changing 
Hq. Ekstein [22] proved that the existence of such operators is a necessary and a sufficient 
condition for two Hamiltonians to be scattering equivalent. 

Any experiment that only measures scattering matrix elements cannot distinguish differ- 
ent scattering equivalent Hamiltonians. If a dynamical representation of the Poincare group 
has a given kinematic subgroup, it is possible to construct unitary scattering equivalences 
[23] that change the kinematic subgroup to any other kinematic subgroup, so any properties 
attributed to a particular kinematic symmetry are not observable. 

III. ASYMPTOTIC STATES 

The probability amplitude for a transition from the prepared initial state to the measured 
final state can be expressed in terms of a scattering operator S. As discussed in the previous 
section, the scattering probability amplitude is the inner product of two solutions of the time- 
evolution equation with past and future asymptotic conditions evaluated at any common 
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time. We can describe neutrino scattering reactions using an ^'-matrix approach. This was 
first done by Cardall [14] within a field-theoretic framework. 

As an illustrative example we consider a reaction where an initial proton and electron 

interact to produce a linear combination of neutrino mass cigcnstatcs which propagate over 
a macroscopic distance until they are absorbed by a second interaction that produces a final 
proton and electron. In the language of second-order perturbation theory this will involve 
two successive reactions, e.g., 

p-\- e n' -\- ^{a) 

n + u -^p' + e'{b). (9) 

This example can be generalized to other cases of interest without affecting the overall 
conclusions, since the role of the intermediate neutrino mass eigenstates is the same. In this 
example, the asymptotic states are given by {p,n,e,p',n',e'} and not the neutrino mass 
eigenstates, which are virtual. 

For the cases of interest, reaction (6) takes place inside a neutrino detector volume, and 
another device records that event via the emerging charged lepton. This is distinct from more 
typical applications in which the S matrix describes one or more events that are distinctly 
separated from any detection equipment. 

The asymptotic free-particle state before the reaction is represented by a localized wave 
packet describing an electron e moving toward the proton p; in addition, there is also a 
free neutron n that is traveling toward the point where it will eventually interact with the 
neutrino mass eigenstates. 

Similarly, the asymptotic free-particle state after the reaction is represented by a localized 
wave packet describing a neutron n' traveling away from the region of the initial interaction, 
and an electron e' and proton p' traveling away from the point where the neutrino mass 
eigenstate was absorbed by the initial neutron n. 

There is a definite probability amplitude for a transition from the initial state containing 
{p,e,n} to the final state containing {p',e',n'}. The unusual feature is that there is not 
a single localized space-time region where the initial and final states overlap. Because we 
only need to construct the scattering states, |^±(i)), whose inner product can be evaluated 
at any common time, the S matrix can be computed in the same manner that is used with 
more traditional asymptotic states. 
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For the example discussed above, the reaction is characterized by two disjoint space-time 
regions locahzed about Xa, where the coherent neutrino mass eigenstate superposition is 
created, and Xb, where it is absorbed. In our example the initial proton p and electron e and 
final neutron n' are localized near Xa at a common time ta, and the initial neutron n and 
final proton p' and electron e' are localized near xi, at a later common time t^. 

In what follows we generalize the expressions to include all lepton fiavors (e, /i, r), labeled 
by Greek indices a, (3. 

The construction of asymptotic states corresponding to this reaction starts with normal- 
izable single-particle states that localize the particles near at time ta or Xf, at time tf, with 
the appropriate expectation value for the initial or final momenta of the observed particles. 
In order to localize these states at the different points x^ and xj we initially localize them 
at the origin and use single particle spatial translations, Tj{—Xc)\4>j), c — a or b, to localize 
the j-th particle near x^ or X;,. The — sign is consistent with 

X, = (0|(X + Xe)|0,) = (0,|T,(Xe)xi;t(x,)|0,.) = (0,. jl^t (_x,)xT,- (-X,) |0) . (10) 

Prom Eq. 4, we have that 

|*±(t)) = ^^±|$±(t)). (11) 

This means that the initial and final scattering states can be transformed to a common time 
by transforming the non-interacting asymptotic states to the common time. To construct 
the non-interacting multi-particle asymptotic states at a common time, all of the particles 
that are near x& at time tb are time-translated to where they would be at time t — ta using 
single-particle time evolution: 

Tj{ta - tb)Tj{-^bM) = Tj{ta - tb, -y^bMj). (12) 

Choosing ta — gives initial and final time-zero asymptotic states of the form 

l*a(0)) = llT{-XaMa)llT{-Xt,)\<Pja). (13) 

i j 

Here the single-particle states \(j)ia) arc initially localized at the origin, and all of the "reac- 
tion geometry" appears in the one-body space-time translation operators, Ti{—Xi). In the 
absence of interactions these non- interacting asymptotic states, |$a(0)), put the initial or 
final reaction products that interact at time t — ta — near x^ and the initial and final re- 
action products that interact at time t — ti, near x^. Because all of the "reaction geometry" 
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is in the structure of |$a(0)), the formal scattering operator, S, is unchanged and can be 
calculated using standard methods of quantum mechanics or quantum field theory. 

In what follows we use the notation Sfi{xb,Xa) for iS'-matrix elements with asymptotic 
states describing neutrinos created at space-time point Xa and absorbed at space-time point 

Xb- 

In this formulation of the scattering problem, the space-time displacement, xi,a '■— Xb — Xa, 
transforms like a four- vector under Lorentz transformations. To show this, note that the 
scattering matrix element is Lorentz invariant: 

l*a) ^ l*a') U{A, 0)\^t) ^ Sf, = S'f,. (14) 

The effect of this transformation on the asymptotic states follows from the Poincare invari- 
ance of the wave operators 

U(A,a)n±^n±Uo(A,a) (15) 

where Uo{A, a) is a product of single-particle unitary representations of the Poincare group, 
Uj{A,a). The Poincare invariance of the wave operators ensures that the representations 
of the Poincare group for the free and interacting systems agree when the particles are 
asymptotically separated. Because Tj{xc) — Uj{I,Xc), the group representation property 
implies 

Uj{A,0)Tj{-x,) = T^{-Ax,)Uj{A,Q). (16) 
It follows from (15) and (16) that 

Sfi{Xb,Xa) = Sfi>{AXa,AXb) (17) 

where the single particle states localized at the origin are replaced by the transformed states 
localized at the origin 

|0;j = C/,(A,O)|0,-«). (18) 

Thus, ill this formalism, dynamical Poincare traiisformatioiis imply the expected transfor- 
mation properties of the parameters that describe the geometry of the interaction region. 
This will be relevant for the oscillation formula, which in this formalism is determined by 
considering how the scattering operator behaves when these parameters are varied. 
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IV. SECOND-ORDER S MATRIX 

In this section we examine the conditions that lead to neutrino oscillations in the 5"- 
matrix approach. We assume that the weak interaction can be treated perturbatively. In 
our example the leading contribution to the scattering matrix is of second order in the weak 
interaction. As noted above, the initial weak interaction at space-time point Xa produces 
a superposition of neutrino (or anti-ncutrino) mass eigenstates, which propagate to the 
space-time point Xf, of the final weak interaction. We assume that the initial particle(s) 
collide or decay near the space-time point Xa producing a neutrino or anti-neutrino and final 
particle(s). The neutrino or anti- neutrino travels and interacts with initial particle(s) near 
space-time point X}, to produce final particles. The superposition of neutrino mass eigenstates 
is never directly observed; it appears only as an intermediate state. The second-order 
calculation can formally be done in the interaction representation, where the Hq includes all 
interactions except for the interactions with the neutrinos. 

For the purpose of illustration, we will assume that the kinematic subgroup of the Poincare 
group is the three-dimensional Euclidean group, generated by spatial translations and rota- 
tions. This is called an "instant form" dynamics. In an instant-form dynamics it is natural 
to label single-particle basis vectors, |p), by the particle's three momentum p, while the 
fourth component satisfies the mass-shell condition = (m^ + p^)^. Spin labels are sup- 
pressed. The weak interaction, V, that couples to the neutrinos creates fiavor eigenstates 
that can be decomposed into mass eigenstates that propagate. In an instant-form dynamics 
the neutrino production interaction is invariant with respect to kinematic translations and 
thus conserves the sum of the single-particle momenta. As mentioned earlier, this conser- 
vation law is related to the freedom to choose among scattering-equivalent interactions and 
does not affect the S matrix. There are many other scattering-equivalent interactions that 
do not have this kinematic symmetry. 

The second-order contribution to the kernel of scattering operator, S, for our example 
reaction is 





-27Ti6^{pa+Pb-p'a-Pb)^ 
^ Ea- Ea' - E,. + iO+ 



(19) 
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where Pa-iPbiP'a-iP'b initial and final four momenta at points a and h. The Roman 

index j labels the mass eigenstate, and the neutrino energy is 



E^.^^m^.-ipl.. (20) 

In addition, the following relations are specific to an interaction with the three-dimensional 
Euclidean group as a kinematic subgroup: 

(Pna> Pv,|V"ia|(Paa, Ppa)") = 5(p^„ + - Paa - Ppa)(Pna> VviWjaWiVaa, VpaV) ■ (21) 

The sum in (19) is over neutrino mass eigenstates. The ()=*= states are electron- proton 
scattering eigenstates. Equation (19) is a generalization of the standard "two-potential" for- 
mula of Gell-Mann and Goldberger [24] where cluster properties have been used to factorize 
the incoming and outgoing scattering states into products of independent scattering states 
associated with the reactions at Xa and X},. 

The neutrinos enter in the matrix elements 

((P^6,P;6) + ll^,^l|Pn6,P.,) (22) 

and 

(23) 

The three-momentum delta function in Eq. (21) implies that in this representation the 
virtual neutrino three-momentum is constrained by the external kinematics: 

Vuj = Paa - Ppa " Pna- (24) 

The flavor of the initial and final states determines the fiavor of the created/detected 
neutrinos. These matrix elements are constructed from an elementary vertex that involves 
the neutrino mixing matrix. 

It is useful to include explicitly the one-body space-time translation operators used in 

the construction of the non-interacting asymptotic states with the kernel of the scattering 
operator, S. Using the four-momentum conservation of the second-order contribution to the 
S operator gives the expression 
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{{P'jlr P'phy\\Vll\\Pnh- P.o)(p;„- P,o HV}nJ|(Pn.- P;,.) ) 



A probability amplitude is obtained by integrating this kernel over initial and final wave 
packets associated with states localized at the origin. The relevant matrix elements have 
the form 

Sfi{Xb,Xa) = 
<P*PfiPpb)<P*pfiPpb)<f>nf{Pna)dp(3bdppfdpnaX 

dPaadPpadPnbfp ai{Paa) 4'pi ( Ppa ) 4'ni ( Pnb ) 
-27ri J (f)*pf{p^b)(j)*pf{ppb)(f)lf{Pna)dp0bdppfdpnaX 

5'{paa + Ppa + Pub ' p'f^b ' P'pb ' Je^^^^^^+^P^-^-^^^--^ X 

na' P'^j ll^'alKPaa' Ppa) ) 

2-^ E -\- E — E' — E 4- 70+ 

■ ^aa ^pa ^na ^I'j ' 

dpaadPpadpnb<f>aiiPaa)<f>piiPpa)<f>ni{Pnb) H ■ (26) 

Equation (26) assumes that the weak interaction can be treated using perturbation theory. 
It was derived in a representation in which the kinematic subgroup is the three-dimensional 
Euclidean group. Equivalent calculations using different kinematic subgroups will have a 
similar form, but will have different "off-shell" neutrino variables. This equation provides 
the starting point for derivations of the oscillation formula. 



V. ESSENTIAL INGREDIENTS FOR OSCILLATIONS 

Equation (26) describes the scattering matrix element for production and absorption of 
a superposition of neutrino mass eigenstates as a second-order perturbation in the weak 
interaction. Oscillations require that this expression can be approximated by a sum of 
non- vanishing neutrino mass combinations multiplied by different phase factors. We now 
examine the essential assumptions that are needed to derive the standard simple oscillation 
formula in the literature from Eq. (26). 



11 



A. (almost) ultrarelativistic 



Experiments focused on neutrino mixing lead to constraints upon the squared differences 
among neutrino masses, with a scale much less than that of the virtual neutrino momenta, 
which make the neutrinos very relativistic. 

An important property of very relativistic neutrinos is that neutrinos with different masses 
propagate with approximately the same velocity, near the speed of light. This ensures that 
the wave packets associated different mass neutrinos continue to overlap over macroscopic 
distances where oscillations can be observed. 

For very relativistic neutrinos the energy of a mass eigenstate with three-momentum 
Pj — \pj\ can be approximated by: 

2 

+ (27) 

In a representation like Eq. (26), where the interactions are chosen to be kinematically 
translationally invariant, the neutrino momenta, pj, are the same for all j. Oscillations are 
sensitive to the difference of these energies for different mass neutrinos. Contributions from 
corrections beyond this expression are very small in the relativistic limit. 

The leading phase differences are thus proportional to (m| — m^) = {nii — mj){mi + rrij). 
This needs to be small enough so variations in the phase are small over the scale of the size 
of the single particle wave packets. 

It is not clear whether neutrinos with larger mass differences would lead to oscillations. 
It has been shown, for example, that a non-relativistic massive neutrino cannot oscillate 
against much lighter partners [4]. 

B. interactions and factorization 

The operator Vja in Eq. (26) connects the leptonic and hadronic spaces, but factors into 
separate components for each space. This is the case whether one uses the contact Fermi 
interaction or a diagram with W exchange. The hadronic contribution will involve vector 
and axial vector current matrix elements that are measurable; they represent a separate 
factor in the S matrix and contribute to the normalization of S, but not the oscillation 
formula. 
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The leptonic contribution for our example is a matrix element of the charge current 
interaction of the form This matrix depends in principle on the masses of the 

neutrino j and the charged lepton f3. However, the leading contribution to the neutrino mass 
comes from the oscillation phase described above, and so we can substitute an approximate 
form for the matrix element by substituting Uj — )■ z/q, where here and henceforth denotes 
a massless neutrino of arbitrary flavor. The lepton matrix element does, however, depend 
upon the mass of the charged lepton. Put another way, the response of the detector to a 
specific charged lepton flavor depends upon the kinematics of this matrix element. 

In the absence of mixing, the weak interaction could be expressed as an operator Va for 
each lepton flavor a. With mixing, the weak interaction couples each lepton flavor a to each 
of the neutrino mass cigcnstatcs j via an operator Vja- 

These points permit us to factor the weak interaction in the following schematic way: 

Vja = UjaVoa, (28) 

where U is a unitary mixing matrix associated with the lepton-neutrino side of the interac- 
tion, and Vqq depends upon the mass of the charged lepton but employs a massless neutrino 
of flavor a. This approximation places all of the neutrino mass and mixing information into 
the matrix U. 

Different representations of the Poincare group for this process will in general involve 
different Va that have different kinematic factors. However, the ability to factorize the 
interaction via Eq. (28) will lead to the same Uja- That is, the interaction could have 
a representation-dependent form, but very small neutrino masses lead to the ability to 
factorize and thereby extract a mixing matrix independent of the representation [25] . 

To the extent that particle kinematics can deviate from the mass or energy shell, con- 
tributions from the interactions can enter. This means that in general, approaches using 
different forms of dynamics will yield different results at the level of perturbation theory, 
although the exact expression to all orders will agree. For the problem at hand, off-shell 
neutrino kinematics corresponds to higher order weak contributions, which are very small in 
magnitude; as noted below, the deviations from mass/energy shell kinematics are also very 
small. 

Henceforth we also omit the ± labels in Eq. 26 and use plane- wave lepton and hadron 
states. 
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C. delta functions and packets 



Equation (26) reflects the use of time-ordered perturbation theory, in which the inter- 
mediate neutrinos are on their mass shells and the intermediate state can have an energy 
different from that of the external initial or flnal state. The energy denominator in Eq. (26) 
can be written as 



E 4- E — E' — E -+- ?n+ V "u. . na -^o' ■ ■ E + E — E' — E 

(29) 

The energies of the initial and final particles that appear in this denominator depend on the 
three-momenta of these particles. When these momenta are integrated over the factors in 
the numerator of Eq. (26), the delta function term only gives a contribution for momenta 
that give an on-shell neutrino energy, which leads to the phase factor exp(— i£',^T'), where 
the time difference T — ti, — ta is the macroscopic time difference used to construct the 
asymptotic states, which is very large on microscopic scales. 

The principal-value term depends upon an integral that involves the phase factor 
exp (—iET) where T is large and E depends on the momenta. It is not hard to show 
that if the remaining part of the integrand is a smooth function of E, then as T becomes 
large the principal value contribution will be exponentially suppressed relative to the delta 
function in Eq. (29). In practice the integrand includes Jacobians involving square roots 
and wave packets of unknown smoothness. Less restrictive assumptions on the smoothness 
of the integrand lead to an algebraic suppression of the principal value term for large T, but 
the contribution is still very small in practice. The general case was carefully examined by 
Grimus and Stockinger [6]. 

The dominance of the energy delta function in Eq. (29) means that all four components 
of the neutrino four-momentum are constrained: 

p'u.-pL-p'pa-pL- (30) 

If we use this result to evaluate Eq. (26) using plane-wave momentum states, then there is 
no oscillation. To see this, we note that oscillations require the superposition of two or more 
neutrino mass eigenstates. Any particular combination of external particle four-momenta 
will yield at most one neutrino four-momentum p,^^. satisfying Eq. (30) [12, 13], and satisfying 
the mass condition pl^ — m/j of only one mass eigenstate. For plane-wave external particle 
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states, there can be no superposition of neutrino mass eigenstates that would set up an 
interference. 

Thus, the oscillation effect necessarily entails momentum distributions, or packets, of the 
external particles. Under such distributions, Eq. (26) will in general give rise to a superpo- 
sition of neutrino mass eigenstates, with each term carrying slightly different four-momenta. 
Since the energy and three-momentum values will vary inside the external wave packets, so 
too the energy and three-momentum components of the neutrino mass eigenstates will differ 
from each other. The most that we can say is that the contributions from different mass 
eigenstates will have quantitatively similar components of their four-momenta as dictated 
by the momentum ranges of the external wave packets. 

It is important to note here that the term packet refers to the momentum distributions 
of the external particles. The neutrino mass eigenstates also have a momentum distribution 
that will depend upon the external packet shapes via the four-momentum condition of Eq. 30, 
but this distribution is not a wave packet in the usual sense of preparation of observable 
states. 

We have seen that external wave packets are essential to obtaining an oscillation formula. 
At the same time, the relevant distributions must be narrow enough that the phase factor in 
Eq. (26) can be factored out of the integrals over these distributions. When the conditions 
for oscillation are satisfied, the phase is approximately stationary. The external distributions 
must therefore be narrow enough to maintain this property. 

A somewhat curious conclusion to this train of reasoning is that, once the phases in the 
oscillation formula are fixed by the kinematics, they depend, to lowest order, only on the 
neutrino eigenstate masses and the average neutrino energy as determined by the external 
particle kinematics. One can therefore make the following replacement: 

where, as noted above, uq denotes a massless neutrino of arbitrary flavor. 
D. a matter of time 

The expression (26) depends upon the space-time separation between the points (x^, to) 
and (x6,t6). These points could be determined by events that can in principle be measured. 
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However, in a typical oscillation experiment (e.g. a reactor and a detector located a few 
kilometers away), the spatial components Xa and x;, (and thereby the displacement L) are 
essentially known (at least to within a scale characteristic of the reactor size), but the time 
difference T = — is not: the detector could record the time of a neutrino detection 
event, but there is no corresponding record of the initial weak interaction event. Thus, the 
detector in principle is sensitive to a range of possible times T. As we shall see, the range of 
values is in practice quite restricted, to the point where it is possible to use any single value 
of T within this range and obtain the usual oscillation formula. 

We now examine the variation of the S matrix with the time T in Eq. (26). Equivalently, 
wc can consider the variation as a function of a velocity v = |L|/T. Since the neutrino mass 
eigenstates are almost ultrarelativistic, we could first assume that v = c = 1. Alternatively, 
we could assume that v corresponds to the velocity of one of the neutrinos: 

V = V, = |. (32) 

These seem to be reasonable assumptions, given the kinematics. For oscillations it should 

not matter which assumption is made. Furthermore, conditions should be such that the 
oscillating phase approximately factors out the integrals in (26). 

To study these conditions consider the interference phase 0i2 from mass eigenstates 1 
and 2: 



012 ^{Pl-P2)-X^ L 



- (IPiI - IP2I 



(33) 



V 

where the time T has an associated velocity v = L/T. If the neutrino masses mj are small 
compared to their energies, then 

^i^|pil + ^; i^2^|P2| + ^, (34) 

2|Pi| 2IP2I 
Now define p, 5p such that |pi,2| — p±5p. The momentum difference 5p has a maximum 
value determined by the momentum distributions of the external particles. We also express 
the velocity parameter v in terms of an average velocity v as v = v + 6v, where 

The phase is then 



012 ~ -f' 



9 2 / 2 2 

m\ 777.2 I l\ / ^1 ^2 



2|pi| 2IP2I V4|Pir 4|p2 



+ (|Pi|-|P2l) T^ + Trrn^-^^' 



(36) 
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Now define p, Sp such that |pi,2| — p^ Sp. If we further assume that Sp <^ p, then 

2 2 

, mi — rrin 



26p6v + 0{6p'). (37) 



2p 

For the specific case that 5v — 0,we are left with 

<f>i2^L^^^^^ + 0{Sp'). (38) 
2p 

In principle, 6v is arbitrary (as is the time), but its contribution to the phase is limited by a 
factor involving 6p. In fact, 6v is further limited by the range of the pi values that determine 
it. For example, if, at one extreme, v — vi, then 



ml — ml ^ {ml — m^Sp 



2p ■ 2f P"' 

that is, it generates a phase contribution that is suppressed by order 5p/p. Thus, we find 
that 5p <^ 1 in order to have interfering neutrino mass eigenstate contributions, and to fix 
the time T to within correction factors that can be neglected. 



VI. RESULT AND DISCUSSION 

Combining all of the simplifying results of the previous section, we find that the space- 
time phase factors and the mixing matrices factor out of the integrals over external wave 
packets, leaving a scattering matrix element of the form 

Sfi{x,,xa) ~ 5^e-*^"''.-("''-"")t/,-«t/;^s/,(/3,«,j) (40) 
j 

where 

Sfi{/3,a,j) = 

-2tt'^ J dpafdpbf(l>lf{Pbf)<l>*af{Paf)S^{Pa+ Pb - P'a - Pb)x 
{p'l3b^Ppb\\^Ol3\\Pnb,Pi^o)^iEa " K' " Ko){Pna^ PuoWoaWiPaa, Ppa)) X 

dPaidpbi(f)bi{Pbi)(l)ai{Pai)- (41) 

Integrating over wave packets satisfying the conditions described above, wc obtain a 
measurement probability for neutrino production via initial flavor a and neutrino absorption 
via final flavor (5: 

prei ^ \Sfi{xb,Xa)\^ ^ ^ fZ/^e-'^^". •("''-^'^)t/,„t/*„e-*^"''^-(^''-^'^)t/fe, = (42) 

jk 
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S^a + 23fJ 5^ U*pU,^UlJJkf5 exp I 



3>k 




) 



(43) 



where 




2 2 



(44) 



Equation (43) is the standard oscillation formula found in the literature. Its simplicity rests 
upon all of the elements described in the previous section. 

We conclude this section with a brief discussion of what is observable among the elements 
of neutrino oscillations. 

In principle, the mixing matrix U is really part of the interaction Hamiltonian, whose 
matrix elements describe the connection of charged leptons (identified by flavor) with neu- 
trino mass eigenstates, since it is the latter (rather than neutrino flavor eigenstates) that 
must form the basis of a the Poincare group with the correct transformation properties. An 
interaction Hamiltonian in general is not observable, because it is always possible to generate 
scattering equivalent Hamiltonians that leave the observables unchanged. One can constrain 
the problem by completely specifying the Hilbert space, but the extraction will depend upon 
this specification. That said, all of this uncertainty is eliminated if the interaction can be 
treated perturbatively, as it can for many electromagnetic and weak processes. One can still 
perform interaction-dependent unitary transformations that lead to different Hamiltonians, 
but the effect on extracting parameters such as mixing angles will always be of higher order. 

A perturbative ^'-matrix calculation that exhibits dependence upon the choice of kine- 
matic symmetry group in relativistic dynamics is incomplete: the full calculation will be 
unitarily equivalent to other full calculations based upon other symmetry groups. We note 
here that our derivation of the oscillation formula makes use of the symmetry of the three- 
dimensional Euclidean group. This leads to delta functions in three-momenta in the in- 
teraction matrix elements. However, the three-momentum delta functions end up being 
paired (modulo very small corrections) with an energy delta function. As noted above, what 
matters in the end is the resulting four-momentum delta function together with external 
wave packets that determine the kinematics of the neutrino mass eigenstates. One could 
also formulate the derivation using front-form dynamics, which exhibit the symmetry of the 
null plane. The interaction matrix elements contain a delta function in front-form three- 
momenta (p±,p''" = p°-|-p^). In turn, it is paired with a denominator in p" — p'^ —p^, which, 
for the geometry in question, is again approximately a p~ delta function. This again leads 
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to a four-momentum delta function, modulated by external wave packets. 

VII. CONCLUSION 

In this paper we used an /S-matrix approach to study the problem of neutrino oscilla- 
tions. This approach is exact and only involves initial and final states that can in principle 
be prepared and detected in laboratory experiments. From this exact formalism we investi- 
gated the features that are combined to arrive at the standard neutrino oscillation formula, 
Eq. (43). We summarize these briefly below: 

1. The perturbative nature of the weak interaction means that the oscillation formula 
can be derived from a second-order expression via the two-potential equation (26). 

2. All three neutrino masses are small. The leading contribution to the oscillation formula 
comes from the first correction to ultrarelativistic kinematics in the space-time phase. 
All remaining dependence upon neutrino mass can be neglected. 

3. Wave packets are essential for oscillations. The spread of external particle four- 
momenta makes it possible for all intermediate mass eigenstates to contribute co- 
herently to the S matrix with slightly differeing four-momenta. 

4. Explicit off-shell effects can be neglected. Oscillations can be computed in a variety of 
approaches that differ by unitary transformations that depend upon the interaction. 
These approaches will all have the same leading contribution to the oscillation formula, 
with all differences (which can be written as off-shell effects) being of higher order in 
the weak interaction. 

5. The weak Hamiltonians in the second-order calculation can be factored into a product 
of terms, one containing only mixing information and the other describing a weak in- 
teraction involving a fictitious massless, flavored neutrino. Corrections to factorization 
are of higher order in the interaction and/or the neutrino mass. 

6. The oscillation formula depends upon the space-time separation between production 
and disappearance of the neutrino mass eigenstates. The spatial separation is known, 
but the time is not measured. Nevertheless, for reasonably small sizes of the external 
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four-momentum distributions (e.g. detector resolution), the range of relevant times is 
restricted to values that produce the standard oscillation formula. 

We conclude that the existence of a simple neutrino oscillation formula that can be derived 
from a variety of theoretical perspectives represents a remarkable convergence of several 
specific physical properties of neutrinos and their interactions. However, to improve on this 
simple formula will in principle bring back all of the interrelated effects that have so far 
been neglected, and "corrected" oscillation formulae may acquire theoretical dependencies 
that the lowest-order expression does not carry. 

The authors would like to express our gratitude to Boris Kayser for helpful correspondence 
related to this manuscript. One of the authors (W.P.) would hke to acknowledge support 
for this work by the U. S. Department of Energy, Office of Nuclear Physics, under contract 
No. DE-FG02-86ER40286. 
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Neutrino oscillations are examined under the broad requirements of Poincare-invariant scattering 
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I. INTRODUCTION 



There is now strong experimental evidence that neutrino flavor eigenstates are mixed via 
a non-diagonal matrix that connects them to neutrino mass eigenstates [1, 2]. The fact that 
space-time propagation is governed by mass rather than flavor eigenstates gives rise to the 
possibility of oscillation between flavor states over space and time. 

There is a large body of literature addressing various aspects of the quantum mechanics 
of neutrino oscillation. Many of the primary issues were set forth by Kayser [3] in 1981. 
Since then, this subject has been examined from a variety of perspectives [4-16]. 

The phenomenon of neutrino oscillations seems simple enough from the perspective of 
introductory quantum mechanics, for which there are countless examples of oscillations in 
two- and three-level systems. There are, however, differences in the neutrino case that tend 
to work against the intuition of the standard examples: the initial conditions for neutrino 
production in reactors or the Sun are generally not controlled; the coherence length over 
which states are virtual can be kilometers rather than subatomic distances, and this affects 
considerations such as what is large or small in a calculation; one of the weak interactions 
that determine the oscillation interval takes place inside the volume of the detector rather 
than far from it; the neutrino kinematics are ultrarelativistic, but not fully so. There is 
also the distinction between what is measured in an experiment and features of a specific 
theoretical approach. These issues have led to a variety of approaches as cited above. Some 
of these approaches have differing perspectives [12, 14], yet all lead to the same simple 
oscillation formula. 

The goal of this paper is to summarize in compact form the most general principles and 
the necessary physical conditions that lead to the standard neutrino oscillation formula. 

First, we provide an approach to neutrino oscillations from the perspective of scattering 
theory in relativistic quantum mechanics. The oscillation formula is a ratio of probabilities 
that can be obtained in turn via the associated scattering probability amplitudes, i.e. the 
S matrix. The virtue of this approach is its generality and its focus on what can actually 
be measured in oscillation experiments. The S'-matrix approach discussed in this paper has 
been examined [17] [18] using quantum field theory, but the conclusions that follow from 
our approach are not limited to a field theoretic treatment. Our goal here is to provide a 
very general framework utilizing the S matrix and the general requirements of relativity and 
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quantum mechanics. 

Second, we discuss the necessary physical conditions that permit us to use this general 
framework to obtain the oscillation formula. In the process we identify what is, or is not, 
essential for oscillations to occur. 

There are now many papers in the literature devoted to neutrino oscillations that touch 
upon many points addressed in this paper. Our purpose is to cast the problem into a more 
general framework and to enumerate the physical conditions in compact form. 

II. POINCARE INVARIANT QUANTUM MECHANICAL SCATTERING 

It is well known [3] that the standard neutrino oscillation formula depends both upon 
relativity (very low mass neutrino mass eigenstates) and quantum mechanics (interference 
leading to oscillation). For oscillation phenomena in which only the lepton sector partic- 
ipates, a very convenient way to obtain oscillations is to employ electroweak field theory. 
For oscillation phenomena involving nucleons or nuclei, a field theory requires extra care 
since the strong interaction is non-perturbative. In the following discussion we show that 
a standard second-order perturbation expression results from basic properties of relativity 
and quantum mechanics, whether or not one uses a field theory. 

We present here an overview of the basic ingredients of quantum mechanical scattering 
that satisfy relativistic invariance. Further details can be found in Ref. [19]. 

A. relativistic invariance 

Relativistic invariance in a quantum theory means that a change of inertial coordinate 
system is a symmetry of the theory. In special relativity, inertial coordinate systems are 
related by space-time translations and proper orthochronous Lorentz transformations. The 
group generated by these transformations is the Poincare group. Proper orthochronous 
Lorentz transformations do not include the discrete transformations associated space reflec- 
tion and time reversal, which are broken by the weak interaction. Wigner proved [20] that 
the existence of a unitary representation, f/(A,a) (where A represents a Lorentz transfor- 
mation and a a space-time shift) of the Poincare group is both a necessary and a sufficient 
condition for a quantum theory to be relativistically invariant. 
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The Hamiltonian H contains interaction-dependent terms, and generates time transla- 
tions. Consistency of the initial value problem requires that interactions must also appear 
in additional Poincare generators. This means that at most a subgroup of the Poincare 
group can be independent of interactions. When such a subgroup exists, it is called a kine- 
matic subgroup. The largest kinematic subgroups were classified by Dirac [21]. They are 
the three-dimensional Euclidean group (instant-form kinematics), the Lorentz group (point- 
form kinematics), and the subgroup that leaves a plane tangent to the light cone invariant 
(front-form kinematics). Time-ordered perturbation methods stemming from a field theory 
implicitly employ an instant-form representation, but the front form and the point form 
have also been frequently employed in the literature. While the existence of a kinematic 
subgroup is a consequence of the choice of representation of the dynamics, the choice of rep- 
resentation does not change the S matrix. It follows that the choice of kinematic symmetries 
of the interactions have no impact on the observables of the theory; this is an example of a 
representation-dependent feature of an interaction that is not experimentally observable. In 
the context of neutrino physics, interactions in an instant-form of the dynamics are invariant 
with respect to kinematic translations and thus conserve the total three momentum. This 
conservation law does not hold in the other two forms of the dynamics, but the interaction 
in the three different forms are related by unitary transformation that do not change the S 
matrix, which means that the kinematic translational invariance of interactions in an instant 
form of the dynamics is not observable. 

Expressions for the time translation of a quantum system may not exhibit manifest 
covariance, but the existence of a unitary representation of the Poincare group implies that 
all physically observable quantities will have the correct transformation properties. 

B. the S matrix 

Scattering theory in a Poincare invariant quantum theory can be formulated using time- 
dependent methods, based on dynamical and asymptotic Hamiltonians, H and Hq. 

Scattering states {"^(t)) can be expressed in terms of a unitary transformation (time 
evolution) in terms of the full interaction: 

\^±it)) = e-^*|vl/±(0)), (1) 
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where the initial condition, |\l/±(0)), is determined by an asymptotic condition. There are 
two natural scattering asymptotic conditions; they require that the state \'^±{t)) approach 
a state of non-interacting particles, in the asymptotic future or past: 

hm |||vl/±(t))-|$±(t))||=0. (2) 

The asymptotic states have the form, 

|<l>±(t)) = ne-^^°*|<l>±(0)). (3) 

where 11 is a mapping from the Hilbert space of scattering asymptotes to the physical Hilbert 
space [22] that includes the internal structure of the asymptotic particles. This mapping is 
needed when one expresses the weak interaction as an operator in the space of particles that 
exist within a composite system (e.g. quarks in nucleons, nucleons in nuclei). In the particle 
case n is a projection operator while in the field-theory case 11 is constructed in Ref. [23]. 

Equations (1-2) imply that the initial condition for the scattering states are related to 
the initial conditions for the non-interacting states by 

|*±(0)) = lim e'^'Ue-'^°'\^±{0)) := 1^±|$±(0)), (4) 

which also serves to define the wave operators Q^-. 

The probability amplitude that a state prepared to become |<l>_(t)) in the asymptotic 
past scatters into a state that becomes |$+(t)) in the asymptotic future is 

{S) = {^4m-{t)) = {^40)1^40)) = (<fH-(o)|^|$_(o)) (5) 

where 

S := QIQ^ (6) 

is the scattering operator. We have exploited the invariance of the probability amplitudes 
under time translation to emphasize that the probability amplitude can be computed using 
states at any common time. We will make use of this to compute scattering amplitudes 
involving neutrino intermediate states. 

The Poincare invariance of the ^-matrix is a consequence of the Poincare invariance of 
the dynamical theory and the intertwining relations of the wave operators 



U{A,a)n± = n±Uo{A,a) 
5 



(7) 



that relate the dynamical representation of the Poincare group on the dynamical Hilbert 
space, with the asymptotic representation on the asymptotic Hilbert space, "Ho 

This formalism can be used with Hamiltonians having any kinematic symmetry; for any 
Hamiltonian with a given kinematic symmetry there exist equivalent Hamiltonian's with 
different kinematic symmetries that give the same scattering operator. 

C. example 

As an illustrative example we consider a reaction where an initial proton and electron 
interact to produce a linear combination of neutrino mass eigenstates which propagate over 
a macroscopic distance until they are absorbed by a second interaction that produces a final 
proton and electron. In the language of second-order perturbation theory this will involve 
two successive reactions, e.g., 

p + e n' + u (a) 

n + v ^p' + e'{b). (8) 

This example can be generalized to other cases of interest without affecting the overall 
conclusions, since the role of the intermediate neutrino mass eigenstates is the same. In this 
example, the asymptotic states are given by {p,n,e,p',n',e'} and not the neutrino mass 
eigenstates, which are virtual. 

For the cases of interest, reaction (6) takes place inside a neutrino detector volume, and 
another device records that event via the emerging charged lepton. This is distinct from more 
typical applications in which the S matrix describes one or more events that are distinctly 
separated from any detection equipment. The reaction (a) takes place in a disjoint spacetime 
region where the neutrino is produced. 

The asymptotic free-particle state before the reaction is represented by a localized wave 
packet describing an electron e moving toward the proton p near the point Xa where the 
neutrino is produced; in addition, there is also a free neutron n that is traveling toward the 
point Xf, where it will eventually interact with the neutrino mass eigenstates. 

Similarly, the asymptotic free-particle state after the reaction is represented by a localized 
wave packet describing a neutron n' traveling away from the region where the neutrino was 
initially produced, and an electron e' and proton p' traveling away from the point where the 
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neutrino mass eigenstate was absorbed by the initial neutron n. 

There is a definite probabihty amphtude for a transition from the initial state containing 
{p,e,n} to the final state containing {p',e',n'}. The unusual feature is that there is not a 
single localized space-time region where the initial and final states overlap. Because we only 
need to construct scattering states, |\l'±(t)), whose inner product can be evaluated at any 
single common time, the S matrix can be computed in the same manner that is used with 
more traditional asymptotic states. 

D. space-time translations 

The scattering matrix can be calculated using standard methods by including the geom- 
etry of the neutrino experiment in the structure of the asymptotic states. 

We illustrate the construction of the asymptotic states for the example discussed above. 
The reaction is characterized by two disjoint space-time regions localized about Xa, where the 
coherent neutrino mass eigenstate superposition is produced, and Xb, where it is absorbed. 
In our example the initial proton p and electron e and final neutron n' are localized near Xa 
at a characteristic time ta, and the initial neutron n and final proton p' and electron e' are 
localized near Xb at a later characteristic time tf,. The role of the characteristic times ta and 
th is discussed further below. 

The asymptotic Hilbert space, "Hq is the tensor product of single-particle Hilbert spaces. 
The construction of asymptotic states corresponding to this reaction starts with normalizable 
single-particle states that localize the particles near Xq at time ta or x^ at time tb with the 
appropriate expectation value for the initial or final momenta of the observed particles. In 
order to localize these states at the different points, x^ and x;,, we initially localize them 
at the origin and then use single-particle space-time translations Uj{I, (x)) = exp{—ipj ■ x), 
where is the four momentum operator of the j^^ asymptotic particle. The initial and final 
asymptotic states have the form 

\<^±{t)) = n?7o(t)e-^^^'^--^''e-^^P^v^i>|<|(0)) (9) 

which places the initial and final particles associated with neutrino production at time t = ta 
near x^ and the initial and final particles associated with neutrino annihilation at time t — tb 
near x;,. The state |$(0)) corresponds to a product of single-particle wave packets with each 
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one centered near the origin with the appropriate expectation value of momentum. In what 
follows we use the notation Sfi{xb,Xa) for the S matrix constructed from the asymptotic 
states (9). The points Xa and Xf, do not correspond to the precise spacetime location of 
each neutrino production or absorption event, but rather, e.g. the location of the center 
of a reactor or a detector. They are, however, physical space-time points and transform as 
four-vectors under Lorentz transformations. The transformation property 



is a consequence of Eq. (7). Significant variations of space-time production/absorption loca- 
tions about a chosen Xa or Xb will lead to a loss of coherence in the interference amplitudes. 
This is especially true of the time components, and is discussed further below. 

E. second-order S matrix 

In this section we examine the conditions that lead to neutrino oscillations in the S'-matrix 
approach. Since the burden of providing a description of the experiment is taken up by the 
construction of the asymptotic states, the scattering operator or wave operators can be 
calculated using conventional methods. Initially we only assume that the weak interaction 
can be treated perturbatively. 

In our example the leading contribution to the scattering matrix is of second order in 
the weak interaction. As noted above, the initial weak interaction at space-time point Xa 
produces a superposition of neutrino (or anti-neutrino) mass eigenstates, which propagate to 
the space-time point Xf, of the final weak interaction. We assume that the initial particle(s) 
collide or decay near the space-time point Xa producing a neutrino or anti-neutrino and 
final particle(s). The neutrino or anti-neutrino travels and interacts with initial particle(s) 
near space-time point Xb to produce final particles. The superposition of neutrino mass 
eigenstates is never directly observed; it appears only as an intermediate state. The second- 
order calculation can formally be performed using time-ordered perturbation theory (the 
Dyson expansion), which can be used in both the quantum mechanical and field theoretic 
formulations of the problem, expressing the S operator or the wave operators, Q±, as time- 
ordered exponentials of integrals of an interaction term V(t) = Vi(t) + V2{t), where Vi{t) 
includes all interactions except for the interactions with the neutrinos, which appear in V2{t). 




(10) 
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The relevant second-order calculation is obtained by formally treating Vi{t) to all orders and 
V2{t) to second order. The result of the non-perturbative part of the calculation gives initial 
and final wave operators associated with the interaction Vi and requires that the intermediate 
states in the perturbative part of the calculation are eigenstates of Hamiltonian associated 
with the interaction Vi. In principle these wave operators involve composite asymptotic 
states (nucleons rather than quarks). 

For the purpose of illustration, we will use a representation in which the kinematic sub- 
group of the Poincare group is the three-dimensional Euclidean group, which means that 
V{t) is kinematically rotationally and translationally invariant. This is called an "instant 
form" dynamics. In this form it is natural to label single-particle basis vectors, |p), by the 
particle's three momentum p, while the fourth component p° satisfies the mass-shell con- 
dition p° = {m? + p^)^. Spin labels are suppressed. The weak interaction, V2, couples to 
the neutrinos that are flavor eigenstates in the absence of mixing, but instead involves lin- 
ear combinations of mass eigenstates that propagate in space and time. In an instant-form 
dynamics the neutrino production interaction is invariant with respect to kinematic trans- 
lations and thus conserves the sum of the single-particle momenta. As mentioned earlier, 
this conservation law follows from our choice of kinematic subgroup, which does not affect 
the S matrix. There are many other interactions that give the same scattering matrix that 
do not have this kinematic symmetry. 

The second-order contribution to the kernel of scattering operator, S, for our example 
reaction is 

(P^fe) Ppfe) Pnal'S'lPaa, Ppa, Pnb) = 

-2Tri5^{pa+Pb-p'a-p'b)y< 

((pW' Ppby\\VjU\P^b, Pu,){p'na^ Pu,\\Vja\\{Paa, Ppa)~) 

where Pa^Pb^p'aiP'b the initial and final four momenta at points a and h. In this equation 
Vja refers to the operator V2, where the Roman index j labels the neutrino mass eigenstate, 
the Greek index a labels the flavor, and the neutrino energy is 



i?., = y^m., - P^. (12) 
In addition, the following relations are specific to an interaction with the three-dimensional 
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Euclidean group as a kinematic subgroup: 

{p'na^Pu,\Vja\{Paa,Ppay) = 5(pL + Pi^, ' Paa - Ppa) {pL, Pu,\\Vja\\{Paa, Ppa)~) ■ (13) 

Equation (11) is a generalization of the standard "two-potential" formula of Gell-Mann and 
Goldberger [24] based on Poincare invariant quantum mechanics or quantum field theory 
where cluster properties have been used to factorize the incoming and outgoing scattering 
states into tensor products of independent scattering states associated with the reactions 
at Xa and Xf,. This is consistent with the assumption that the interaction, Vi, which was 
treated to all orders, is short ranged. 

The neutrinos enter in the matrix elements 

{ip'^b.PU)^\\V,^p\\Pnb,Pu,) (14) 

and 

(15) 

The three-momentum delta function in Eq. (13) implies that in this representation the 
virtual neutrino three-momentum is constrained by the external kinematics: 

Pu, = Paa - Ppa " Pna- (16) 

The flavor of the initial and final states determines the flavor of the produced/detected 
neutrinos. These matrix elements, which involve initial and final hadrons, are linear in an 
elementary vertex that involves the neutrino mixing matrix. 

It is useful to include explicitly the one-body space-time translation operators used in 
the construction of the non-interacting asymptotic states with the kernel of the scattering 
operator, S. Using the four-momentum conservation of the second-order contribution to the 
S operator gives the expression 

e^(P;,.+P;.-P..)-(-'.--)(p^^, P;„ pLI^IP^a, Ppa, Put) = 

-2m6\p^a + Ppa + Pnb - P'pt - pU - X 

-ir^ {jp'lSb^ P'pbyWVlfsWPnb, Pt.,)(Pna^ Pu,\\Vja\\{Paa, Ppa)~) 
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A probability amplitude is obtained by integrating this kernel over initial and final wave 
packets associated with states localized at the origin. The relevant matrix elements have 
the form 

/ (P/3f>) 0p/ (Ppfe) 0n / {Pna)dpiSbdppfdpna X 

e'^^'^^^^'^^-^-^^^'^"-''^^ {p'^„ pL\S\Paa, Ppa, Put) X 

dpaadppadpnb4> 

-2lTi j (p*f}f{ppb)(i)*pf{Ppb)KfiPna)dpf}bdppfdpnaX 

S\paa + Ppa + Pnb ' p'pb " P'pb " p'^y^^'^^^^'v^-^-^^'^^^-^^^ X 
v;^ ((P^b, Ppfe)^ II II Pn6, Pu^ ) {p'na^ Pv, \\Vj 

a II vPaa! Ppa) / 

^ E -I- E — E' — E + ?n+ 

J -^aa \ -^pa -^na ^ 

dpaadppadpnb'^i ai ( Poa ) 0pi ( Ppa ^'Pni (Pnb) + ■■■■ (18) 

Equation (18) assumes that the weak interaction can be treated using perturbation theory. 

This basic result is valid in both the field theoretic and quantum mechanical formulations. 
It was derived in a representation in which the kinematic subgroup is the three-dimensional 
Euclidean group. Equivalent calculations using different kinematic subgroups will have a 
similar form, but will have different "off-shell" neutrino variables. This equation provides 
a universal starting point for derivations of the oscillation formula. In what follows we 
discuss what approximations to the general formula (18) are needed to derive the standard 
oscillation formula. 



III. ESSENTIAL INGREDIENTS FOR OSCILLATIONS 

Equation (18) describes the scattering matrix element for production and absorption of 
a superposition of neutrino mass eigenstates as a second-order perturbation in the weak 
interaction. A process with initial flavor a and final flavor /3 has a probability of observa- 
tion proportional to |S'/3ap, and the quantities of relevance are ratios of these probabilities. 
Oscillations require that this expression in Eq. (18) can be approximated by a sum of non- 
vanishing neutrino mass combinations multiplied by different phase factors. We now examine 
the essential assumptions that are needed to derive the standard simple oscillation formula 
in the literature from Eq. (18). 
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A. (almost) ultrarelativistic 



Experiments focused on neutrino mixing lead to constraints upon the squared differences 
among neutrino masses, with a scale much less than that of the virtual neutrino momenta, 
which make the neutrinos very relativistic. 

An important property of very relativistic neutrinos is that neutrinos with different masses 
propagate with approximately the same velocity, near the speed of light. This ensures that 
the wave packets associated different mass neutrinos continue to overlap over macroscopic 
distances where oscillations can be observed. 

For very relativistic neutrinos the energy of a mass eigenstate with three-momentum 
Pj = \pj\ can be approximated by: 

2 

171 ■ 

+ (19) 

In a representation like Eq. (18), where the interactions are chosen to be kinematically 
translationally invariant, the neutrino momenta, pj, are the same for all j. Oscillations are 
sensitive to the difference of these energies for different mass neutrinos. Contributions from 
corrections beyond this expression are very small in the relativistic limit. 

B. interactions and factorization 

The operator Vja in Eq. (18) connects the leptonic and hadronic spaces, but factors into 
separate components for each space. This is the case whether one uses the contact Fermi 
interaction or a diagram with W exchange. The hadronic contribution will involve vector 
and axial vector current matrix elements that are measurable; they represent a separate 
factor in the S matrix and contribute to the normalization of S, but not the oscillation 
formula. 

The leptonic contribution for our example is a matrix element of the charge current 
interaction of the form (/SH^y-All'^j)- This matrix depends in principle on the masses of the 
neutrino j and the charged lepton /3. However, the leading contribution to the neutrino mass 
comes from the oscillation phase described above, and so we can substitute an approximate 
form for the matrix element by substituting Uj — z/q, where z/q here and henceforth denotes 
a massless neutrino of arbitrary flavor. The lepton matrix element does, however, depend 
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upon the mass of the charged lepton. Put another way, the response of the detector to a 
specific charged lepton flavor depends upon the kinematics of this matrix element. 

In the absence of mixing, the weak interaction could be expressed as an operator for 
each lepton flavor a. With mixing, the weak interaction couples each lepton flavor a to each 
of the neutrino mass eigenstates j via an operator Vja- 

These points permit us to approximately factor the weak interaction in the following way: 

Vjo. = UjaVoa, (20) 

where U is a unitary mixing matrix associated with the lepton-neutrino side of the interac- 
tion, and Voa depends upon the mass of the charged lepton, but employs a massless neutrino 
of flavor a. This approximation places all of the neutrino mass and mixing information into 
the matrix U. 

Different representations of the Poincare group for this process will in general involve 
different Va that have different kinematic factors. However, the ability to factorize the 
interaction via Eq. (20) will lead to the same Uja- That is, the interaction could have 
a representation-dependent form, but very small neutrino masses lead to the ability to 
factorize and thereby extract a mixing matrix independent of the representation [17]. 

To the extent that particle kinematics can deviate from the mass or energy shell, con- 
tributions from the interactions can enter. This means that in general, approaches using 
different forms of dynamics will yield different results at the level of perturbation theory, 
although the exact expression to all orders will agree. For the problem at hand, off-shell 
neutrino kinematics corresponds to higher order weak contributions, which are very small in 
magnitude; as noted below, the deviations from mass/energy shell kinematics are also very 
small. 

Put another way, in time-ordered perturbation theory, full relativistic invariance is not 
satisfied to finite order. Lorentz transformations involve unitary operators that depend upon 
the interaction, and this in principle mixes the orders of perturbation. Since the interaction 
in question is weak, the violation of relativistic invariance involves higher order perturbative 
corrections that are very small. 

Since the final-state interactions between the electron and proton do not affect the oscil- 
lation formula, which involve ratios of squares of the S matrix, we also omit the ± labels in 
Eq. (18) and use plane- wave lepton and hadron states. 

13 



C. delta functions and packets 



Equation (18) reflects the use of time-ordered perturbation theory, in which the inter- 
mediate neutrinos are on their mass shells and the intermediate state can have an energy 
different from that of the external initial or final state. The energy denominator in Eq. (18) 
can be written as 

E +E -E' -E +zO+ " + - ^'na - ^-.) + +E ^-E' -E ■ 

(21) 

The energies of the initial and final particles that appear in this denominator depend on the 
three-momenta of these particles. When these momenta are integrated over the factors in 
the numerator of Eq. (18), the delta function term only gives a contribution for momenta 
that give an on-shell neutrino energy, which leads to the phase factor exp{—iEiyT), where 
the time difference T = tf, — is the macroscopic time difference used to construct the 
asymptotic states, which is very large on microscopic scales. 

The principal-value term depends upon an integral that involves the phase factor 
exp {—iET) where T is large and E depends on the momenta. It is not hard to show 
that if the remaining part of the integrand is a smooth function of E, then as T becomes 
large the principal value contribution will be exponentially suppressed relative to the delta 
function in Eq. (21). In practice the integrand includes Jacobians involving square roots 
and wave packets of unknown smoothness. Less restrictive assumptions on the smoothness 
of the integrand lead to an algebraic suppression of the principal value term for large T, but 
the contribution is still very small in practice. The general case was carefully examined by 
Grimus and Stockinger [6]. 

The dominance of the energy delta function in Eq. (21) means that all four components 
of the neutrino four-momentum are constrained: 

P'u, = pL - P'pa - pL- (22) 

If one uses plane- wave momentum states to evaluate Eq. (18), then the presence of these 
delta functions precludes an oscillation. Observables such as a cross section or a detection 
rate (for cases where this is no definable incident beam of initial particles) are proportional 
to the square of the S matrix. Oscillations will depend upon the interference between 
matrix elements of S and S*^ involving different neutrino mass eigenvalues, and in turn these 
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matrix elements contain different four-momentum delta functions that cannot be satisfied 
simultaneously, unless the two mass eigenvalues are the same. The result is an observable 
that is the sum of squares of S'-matrix elements for each mass eigenvalue, with no interference 
terms and therefore no oscillation contributions. Oscillations can only appear when the 
external particle momentum plane wave states are replaced by wave packets that then permit 
the product of delta functions to be satisfied by slightly different four-momenta and different 
neutrino mass eigenvalues. 

Thus, the oscillation effect necessarily entails momentum distributions, or packets, of the 
external particles. Under such distributions, Eq. (18) will in general give rise to a superpo- 
sition of neutrino mass eigenstates, with each term carrying slightly different four-momenta. 
Since the energy and three-momentum values will vary inside the external wave packets, so 
too the energy and three-momentum components of the neutrino mass eigenstates will differ 
from each other. The most that we can say is that the contributions from different mass 
eigenstates will have quantitatively similar components of their four-momenta as dictated 
by the momentum ranges of the external wave packets. 

It is important to note here that the term packet refers to the momentum distributions 
of the external particles. The neutrino mass eigenstates also have a momentum distribution 
that will depend upon the external packet shapes via the four-momentum condition of 
Eq. (22), but this distribution is not a wave packet in the usual sense of preparation of 
observable states. It also depends on assumptions that have no impact on the final S'- 
matrix. 

We have seen that external wave packets are essential to obtaining an oscillation formula. 
At the same time, the relevant distributions must be narrow enough that the phase factor in 
Eq. (18) can be factored out of the integrals over these distributions. When the conditions 
for oscillation are satisfied, the phase is approximately stationary. The external distributions 
must therefore be narrow enough to maintain this property. 

A somewhat curious conclusion to this train of reasoning is that, once the phases in the 
oscillation formula are fixed by the kinematics, they depend, to lowest order, only on the 
neutrino eigenstate masses and the average neutrino energy as determined by the external 
particle kinematics. One can therefore make the following replacement: 

SiEaa + Epa - -E.^^ - E^.) -)■ 5{Eaa + Epa - E'^^ - E^^^), (23) 
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where, as noted above, uq denotes a massless neutrino of arbitrary flavor. 
The oscillations are due to the dependence of the phase e*'-^'"''''^?'''"^"''^'^^'' 



in Eq. (18) on 



the neutrino masses. In the limit that the neutrinos propagate on shell the linear combination 
of four momenta in the exponent become four momentum of the propagating neutrino. As 
will be seen below, the relevant four momenta are those of neutrino mass eigenstates, so the 
phase angle in Eq. (18) becomes 



for a mass rrij. 

D. a matter of time 

The expression (18) depends upon the space-time separation between the points (x^, to) 
and (x6,tf,). These points could be determined by events that can in principle be measured. 
However, in a typical oscillation experiment (e.g. a reactor and a detector located a few 
kilometers away), the spatial components and Xf, (and thereby the displacement L) are 
essentially known (at least to within a scale characteristic of the reactor size), but the time 
difference T = th — ta is not: the detector could record the time of a neutrino detection 
event, but there is no corresponding record of the initial weak interaction event. Thus, the 
detector in principle is sensitive to a range of possible times T. As we shall see, the range of 
values is in practice quite restricted, to the point where it is possible to use any single value 
of T within this range and obtain the usual oscillation formula. 

We now examine the variation of the S matrix with the time T in Eq. (18). Equivalently, 
we can consider the variation as a function of a velocity v = |L|/T. Since the neutrino mass 
eigenstates are almost ultrarelativistic, we could first assume that v = c = 1. Alternatively, 
we could assume that v corresponds to the velocity of one of the neutrinos: 



These seem to be reasonable assumptions, given the kinematics. For oscillations it should 
not matter which assumption is made. Furthermore, conditions should be such that the 
oscillating phase approximately factors out the integrals in (18). 

The square of the scattering matrix S^a will involve interference phases — (f)j — 0^, 
where (pj is a contributing phase as defined in Eq. (24). Consider the interference phase (pu 




(24) 




(25) 
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from mass eigenstates 1 and 2: 

(Pi -P2) -x 



>\2 



El — E2 



(IPi 



IP2I 



(26) 



where the time T has an associated velocity v = L/T. If the neutrino masses rrii are small 
compared to their energies, then 



El ^ |pi| 



m^ 



E,. 



IP2I 



rrin 



(27) 



2|pii' ^ ■ 2|p2r 

Now define p, 6p such that |pi,2| = p^Sp. The momentum difference 6p has a maximum 
value determined by the momentum distributions of the external particles. We also express 
the velocity parameter v in terms of an average velocity v a.s v = v + 6v, where 



V = - ( — + ^ 

2 \Ei E2 



rrin 



4|pi 



4|P2| 



The phase is then 



n2 



(IpiI - IP2I) 



+ 



6v 



_2|pi| 2IP2I ■ V4|pip ■ 4|p2|2 

Now define p, 6p such that |pi,2| = p± Sp. If we further assume that 6p ^ p, then 



-'12 



2 2 

^ 

2p 



-2(5p(5w + 0(5p2). 



For the specific case that 6v = 0,we are left with 



^12 



L^^i^ + 0(5/). 



2p 



(28) 



(29) 



(30) 



(31) 



In principle, 6v is arbitrary (as is the time), but its contribution to the phase is limited by a 
factor involving 6p. In fact, 6v is further limited by the range of the pi values that determine 
it. For example, if, at one extreme, v = vi, then 



'12 



[mi 



ml)Sp 



+ 0{6p'), 



(32) 



2p 2p2 

that is, it generates a phase contribution that is suppressed by order 6p/p. Thus, we find 
that 6p <^ 1 in order to have interfering neutrino mass eigenstate contributions, and to fix 
the time T to within correction factors that can be neglected. 

This result is effective for small neutrino masses. It is not clear whether neutrinos with 
larger mass differences would lead to oscillations. It has been shown, for example, that a 
non-relativistic massive neutrino cannot oscillate against much lighter partners [4]. 
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IV. RESULT AND DISCUSSION 



Combining all of the simplifying results of the previous section, we find that the space- 
time phase factors and the mixing matrices factor out of the integrals over external wave 
packets, leaving a scattering matrix element of the form 

j 

where 

Sfi{l3,a,j) = 

-2tt'^ j dpafdpbf(Plf{Pbf)<P*af{Paf)S'^{Pa + Pb - p'a ' P'b) >< 

)5{Ea - Ea' - -E^o)(Pna5P!^ollK)Q||(Paa,Ppa))x 
dpaidVbi(t)bi{Pbi)(t)ai{Pai)- (34) 

Integrating over wave packets satisfying the conditions described above, we obtain a 
measurement probability for neutrino production via initial flavor a and neutrino absorption 
via final flavor /3: 

+ 2M 5^ U*pU,^UlJJup exp (^^^) (36) 

j>k 

where 

Am^^ = - ml. (37) 

Equation (36) is the standard oscillation formula found in the literature. Its simplicity rests 
upon all of the elements described in the previous section. 

We conclude this section with a brief discussion of what is observable among the elements 
of neutrino oscillations. 

In principle, the mixing matrix U is really part of the interaction Hamiltonian, whose 
matrix elements describe the connection of charged leptons (identified by flavor) with neu- 
trino mass eigenstates, since it is the latter (rather than neutrino flavor eigenstates) that 
must form the basis of a the Poincare group with the correct transformation properties. 
An interaction Hamiltonian in general is not observable, because it is always possible to 
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generate scattering equivalent Hamiltonians that leave the observables unchanged. That 
said, all of this uncertainty is eliminated if the interaction can be treated perturbatively, 
as it can for many electromagnetic and weak processes. One can still perform interaction- 
dependent unitary transformations that lead to different Hamiltonians, but the effect on 
extracting parameters such as mixing angles will always be of higher order in the neutrino 
mass differences. 

A perturbative S'-matrix calculation that exhibits dependence upon the choice of kine- 
matic symmetry group in relativistic dynamics is incomplete: the full calculation will be 
unitarily equivalent to other full calculations based upon other kinematic symmetry groups. 
We note here that our derivation of the oscillation formula makes use of the symmetry of the 
three-dimensional Euclidean group. This leads to delta functions in three-momenta in the 
interaction matrix elements. However, the three-momentum delta functions end up being 
paired (modulo very small corrections) with an energy delta function. As noted above, what 
matters in the end is the resulting four-momentum delta function together with external 
wave packets that determine the kinematics of the neutrino mass eigenstates. One could 
also formulate the derivation using front-form dynamics, which exhibit the symmetry of the 
null plane. The interaction matrix elements contain a delta function in front-form three- 
momenta {p±,p'^ = + p^). In turn, it is paired with a denominator in p~ = p^ —p^, which, 
for the geometry in question, is again approximately a p~ delta function. This again leads 
to a four-momentum delta function, modulated by external wave packets. 

V. CONCLUSION 

In this paper we used an S'-matrix approach to study the problem of neutrino oscillations. 
This approach is exact and only involves initial and final states that can in principle be 
prepared and detected in laboratory experiments. We have used very general principles of 
relativity and quantum mechanics to obtain an S matrix and a general expression in second- 
order perturbation theory in the weak interaction. Here we summarize briefly the conditions 
that must be satisfied for the second-order S matrix to yield the standard neutrino oscillation 
formula of Eq. (36). 

1. The perturbative nature of the weak interaction means that the oscillation formula 
can be derived from a second-order expression via the two-potential equation (18). 
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2. All three neutrino masses are small. The leading contribution to the oscillation formula 
comes from the first correction to ultrarelativistic kinematics in the space-time phase. 
All remaining dependence upon neutrino mass can be neglected. 

3. Wave packets are essential for oscillations. The spread of external particle four- 
momenta makes it possible for all intermediate mass eigenstates to contribute co- 
herently to the S matrix with slightly differing four-momenta. 

4. Explicit off-shell effects can be neglected. Oscillations can be computed in a variety of 
approaches that differ by unitary transformations that depend upon the interaction. 
These approaches will all have the same leading contribution to the oscillation formula, 
with all differences (which can be written as off-shell effects) being of higher order in 
the weak interaction. 

5. The weak Hamiltonians in the second-order calculation can be factored into a product 
of terms, one containing only mixing information and the other describing a weak in- 
teraction involving a fictitious massless, flavored neutrino. Corrections to factorization 
are of higher order in the interaction and/or the neutrino mass. 

6. The oscillation formula depends upon the space-time separation between production 
and disappearance of the neutrino mass eigenstates. The spatial separation is known, 
but the time is not measured. Nevertheless, for reasonably small sizes of the external 
four-momentum distributions (e.g. detector resolution), the range of relevant times is 
restricted to values that produce the standard oscillation formula. 

We conclude that the existence of a simple neutrino oscillation formula that can be 
derived from a variety of theoretical perspectives represents a remarkable convergence of 
several specific physical properties of neutrinos and their interactions. However, to improve 
on this simple formula will in principle bring back all of the interrelated effects that have so 
far been neglected, and "corrected" oscillation formula may acquire theoretical dependencies 
that the lowest-order expression does not carry. 
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